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PACS 71.70.Ej - Spin-orbit coupling, Zeeman and Stark splitting, Jahn- Teller effect 

PACS 72 . 25 . -b - Spin polarized transport 

PACS 11 . 15. Ex - Spontaneous breaking of gauge symmetries 

PACS 73 . 43 . -f - Quantum Hall effects 

Abstract - The Hamiltonian describing 2D electron gas, in a spin-orbit active medium, can be 
cast into a consistent non-Abelian gauge field theory leading to a proper definition of the spin 
current. The generally advocated gauge symmetric version of the theory results in current densities 
that are gauge covariant, a fact that poses severe concerns on their physical nature. We show 
that in fact the problem demands gauge fixing, leaving no room to ambiguity in the definition of 
physical spin currents. Gauge fixing also allows for polarized edge excitations not present in the 
gauge symmetric case. The scenario here is analogous to that of superconductivity gauge theory. 
We develop a variational formulation that accounts for the constraints between U{1) physical 
fields and SU (2) gauge fields and show that gauge fixing renders a physical matter and radiation 
currents and derive the particular consequences for the Rashba SO interaction. 



Introduction. — In condensed matter physics appli- 
cations, the spin current is often considered an ill-defined 
notion [T], or at least a non conserved quantity since an- 
gular momentum may be transferred via Spin- Orbit inter- 
actions (SOI) to the lattice. Obviously, the total angu- 
lar momentum is conserved so it is illuminating to follow 
the route proposed by Yang and Mills in their description 
of the consequences of the isotopic spin conservation [2]: 
The conservation of isotopic spin points to the existence of 
a fundamental invariance law similar to the conservation 
of electric charge. In the latter case, the electric charge 
serves as a source of electromagnetic field. An important 
concept in this case is gauge invariance which is closely 
connected with (1) the equation of motion of the electro- 
magnetic field, (2) the existence of a current density, and 
(3) the possible interactions between a charged field and 
the electromagnetic field. In the present context, we es- 
sentially start from point (3), since we build a consistent 
gauge theory from the SOI (see Refs. [3HI2]) and we de- 
duce point (2), i.e. the definition of a conserved quantity, 
the dissipationless spin current density. A problem then 
arises, due to the non-Abelian character of the underly- 



ing gauge theory. The current density (color current in 
the context of high energy physics, but here an angular 
momentum current density) consists of two parts: a mat- 
ter contribution plus the so-called radiation contribution, 
the existence of which is linked to the fact that the non- 
Abelian gauge fields carry the corresponding color charge 
(this is also the origin of the non linear character of the 
equations of motion satisfied by the gauge fields, point (1) 
in Yang and Mills' quotation). The matter contribution 
here corresponds to the spin current and the radiation 
contribution is associated to the angular momentum car- 
ried by the lattice degrees of freedom. The difficulty is 
that the conserved current density has atrocious transfor- 
mation properties under a change of gauge (Ramond [13] , 
explicitely given in Ref . [14j ) , and in the gauge symmetric 
formulation [3 [11], neither the total current density nor 
the two contributions separately (matter or radiation) is 
gauge invariant. None of these densities is then acceptable 
as a physical quantity. 

In this paper we draw attention to the fact that both the 
Pauli equation and the gauge formulations of the Rashba 
and 2D Dresselhaus Hamiltonian cannot be written cor- 
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rectly in a generally gauge invariant fashion [HUH]- We 
argue that this makes a very strong physical statement, 
since as discussed above, a gauge invariant Hamiltonian 
would not allow a locally gauge invariant spin current or 
spin polarization. We show that, in fact, the correct for- 
mulation contains and demands gauge symmetry breaking 
that allows for non ambiguous matter spin current. This 
is also related to the existence of a gauge potential as a 
physical field in the same way as in superconductivity, re- 
quired by the very identification made in where it is 
directly related to a physical U{1) electric field. Note that 
from a semantic point of view, we will follow Sonin's use 
of the term superfluidity, superfluidity means only a possi- 
bility to transport a physical quantity ( mass, charge, spin, 
etc.) without dissipation [T] . 

Yang-Mills formulation. — We first recall the es- 
sential steps of the construction of the gauge field theory 
associated to SOI. Let us start by considering the Pauli 
Hamiltonian to order v/c (in the International System of 
units) 



H 



(P - eA) 



V 



eh 



2m 2m 
ehcr • E X (p — eA) 



B 



,h' 



Am^c^ 
rV -E, 



3?i^cr • V X E 



(1) 



where the first and second terms are the kinetic energy, in- 
cluding the minimal coupling to the electromagnetic field 
(e = — |e| is the electron charge) and the substrate po- 
tential denoted by V , that can be assumed periodic. The 
third term is the Zeeman interaction where B is the mag- 
netic field and cr is the vector of Pauli matrices cr = ct^Uq. 
The fourth and fifth terms are the spin-orbit interaction 
and the last one is the Darwin term (E is the electric field 
and the bold font is used to denote vectors in ordinary 
space). 

In most of the applications considered, the rotor of the 
electric field is absent (TS], however, as it is discussed be- 
low, our main results would not depend on this assump- 
tion. In order to suggest an SU{2) x U{1) form we can 
rewrite the Hamiltonian, following Jin, Li and Zhang [7] 
as 



H 



eA 



9 

'Sm 



(2) 



One can identify the third term in the square 
as a new SU{2) connection defined by gW^T"- = 



-{eh/2mc 



or 



{e|2m^)ta3^EJ [IS] 



where the r" (half the Pauli matrices) are the symmetry 
generators for SU{2) obeying the commutation relation 
[t°,t''] = iCabcT'^: the symbol Cabc being the totally anti- 
symmetric tensor, and g = h. The relation between the 



spin operator and the corresponding generators is gr" = 
s" with spin one-half. The fourth term in equation © 
represents the Zeeman energy and defines the time com- 
ponent of the SU{2) connection, gcW^r"" = — (e/m)s • B 
(or Vl^Q = —{e/mc)B°'). Following the conventional use for 
quadri-vectors we set = {Wq, W°^) and = {Aq, Ai) 
with Aq = V/ce in equation ([T]). The Latin indices from 
the middle of the alphabet, ijk, denote the spatial coor- 
dinates and run over three values, Greek indices play the 
same role including also as the time component, and 
the letters a,b,c denote the internal space indices for the 
generators. We always use the internal space indices as 
superindices, except in the antisymmetric tensor. 

This formulation differs crucially from that of ref. [l][7l 
[TT| in the second term written here as a function of the 
SU(2) connection to evidence gauge symmetry breaking 
(GSB) in this Hamiltonian. The origin of this term is 
found in the fact that the SOI is linear in momentum. 

The Lagrangian generating the correct Schrodinger 
equation is thus 



C = 



-V^ W^^W^f + hcW^T" + ecAo ) V 



(3) 



where -0 is a Pauh spinor, GJJj, = d^jW^ — dyW^ — 
eabcW^W^ and F^, = d^A, - d,A^ are the SU{2) and 
U{1) field tensors, respectively. The free non-Abelian 
field contribution —jhcG'^^G"'^'^ corresponds to the spin- 
orbit correction to the energy in a rotationally invariant 
system The covariant derivatives are defined by 

Vi ~ di — ^Ai — iW^T°- for the space derivative and 
T>t = dt + ^Aq + icWqt'^ for the corresponding time 
derivative. The first two lines of equation ^ can be re- 
ferred to as the matter part £mat. (since they contain the 
spinor ip) and the last line as the pure field (or radiation) 
part £fioid- 

Let us illustrate our purpose with the case of a quasi 
two-dimensional non interacting electron gas subject to 
both Rashba and Dressclhaus spin orbit interactions [TT] . 
The components of the non-Abelian gauge field arc in this 
case 



2m 



/3T^)y. 



(4) 



It is possible to measure and control the Rashba pa- 
rameter a using gate voltages, e.g. in two dimensional 
GaAs/AlGaAs electron gas [HHST], a = b{E), where {E) 
is the expectation value of the electric field at the 2DEG, 
and b depends on the inverse of both the effective mass and 
the material gap [22] . The simple form for the Dresselhaus 
term (with parameter /?) is obtained with the assumption 
of strong confinement at the length scale d of the electron 
gas, hp ^ ir/d [53] where kp is the Fermi wa vector of 
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the in plane electrons, reducing cubic contributions to the 
dominant linear terms in equation 

Under a global change of gauge, the gauge field dU 
would rotate in spin space, i.e. the physical situation 
would be different, with a SO interaction resulting from a 
rotated electric field. Such a gauge transformation would 
change the effective magnetic field around which the spin 
precesses, k x E, and change the workings of SO based de- 
vices. The very values of the SOI Rashba and Dresselhaus 
would be changed one into another by such a transforma- 
tion, while experimentally they have definite values. This 
simple physical argument shows that a mechanism should 
exist to prevent from such gauge transformations, and this 
is the subject of the next section. 

Gauge fixing. — It is clear that the Lagrangian den- 
sity ([3]) is not invariant under local S'[/(2)— gauge trans- 
formations U = e^'^ " 



ip —i' Uip, 



(5) 
(6) 

(7) 



The presence of the quadratic, in the gauge field, term 
restricts the gauge invariance to a smaller set of transfor- 
mations. Let us focus the attention to this term. 



(8) 



8m 

Under an infinitesimal local gauge trans- 
formation parameterised by a", this con- 
tributes an additional term in the Langrangian, 
I^V^ {W^d'a" + (aja")t¥") V + 0(a2) which can 
be written as a pure divergence, j^di (M^"a"VV) ~ 

^ ((5iTV")a°i/'V + W^^a^S,; (^AV)) at the condition 
that the gauge transformations are restricted to 



d,W'" = and 9, = 0. 



(9) 



The first condition is an SU{2) version of the Coulomb 
gauge, leading here to V x E = which is consis- 
tent with our choice of simplification of the SO inter- 
action. In a more general situation, one would have 



-(e/mc^)[V X E -I- dt'B]i, = 0, also consis- 



satisfies such a con- 



tent with the Maxwell- Faraday law |8j. The second con- 
straint in equation ([9]) is the condition for fixing the 
norm of the wave function. It is worth noting that a 

spinor of the form i' ^ ^e»e(x) 

ditions and leads to a superfluid matter spin current (de- 
fined later) J° = |^(V0)n(5°''. In the partition function 
Z = J 'D[^p^'D[^p]e~ i' -f this term contributes a high 
weight in the kinetic energy and in the mean field approx- 
imation, the excitations with fluctuating phase can be ne- 
glected, leaving zero bulk matter spin current, a situation 
similar to superconductivity. After the gauge transfor- 
mation, the remaining pure divergence contributes to the 



action as a surface term, leaving the action covariantly 
transformed 



S 



d*xC + ^ J dt j da^W'tt^ipH 



S + surface terms. 



(10) 



In a consistent gauge formulation, the surface term being 
gauge dependent should vanish, leading to boundary con- 
ditions at the edge of a finite system. This would imply 
that the non-Abelian gauge vector has no component per- 
pendicular to the surface, W°|edge'n = 0, Va. This condi- 
tion induces a similar constraint for the "non-Abelian elec- 
tric field" E° to be defined in the following section, and as 
this field will be seen to obey Maxwell-like equations, the 
discontinuity of its normal component across a boundary 
is associated to the existence of a spin polarization sitting 
at the boundary. Vanishing normal components thus do 
not allow for spin accumulation at the edges of a finite sys- 
tem, in contradiction with experimental results \27\ . From 
this argument, one has to conclude that there should be 
no gauge freedom, hence no surface contribution to the 
action. This is obviously corroborated by the fact that 
the non-Abelian gauge fields and Wq are defined in 
terms of the physical U{1) electric and magnetic fields. 

The situation is similar to the case of superconduc- 
tivity, for which there is also a mass term that gen- 
erates, through gauge transformations, gauge-dependent 
bulk and surface contributions. Killing the bulk contribu- 
tion necessitates imposing the Coulomb gauge condition 
V • A = and a condition on the superconducting wave- 
function V|(pp = 0. There still remains a surface term 
which disappears under restricted gauge transformations 
at the condition that A be tangent to the boundary of 
the system. This would then imply vanishing tangential 
components for the magnetic field, hence no surface su- 
perconducting current. This is obviously a wrong conclu- 
sion which can only be discarded if from the beginning one 
would not allow any gauge freedom. This mechanism leads 
to the phase locking of the wave function, and renders the 
gauge vector field physical, since it is automatically pro- 
portional to the superfluid current. 

Returning to the spin-orbit interaction, in order to ob- 
tain the equations of motion for the gauge flelds, we have 
to take into account another peculiarity of the present 
problem, i.e. the fact that the U{1) and the SU{2) gauge 
fields are not independent. The variational formulation 
has to include constraints via Lagrange multipliers which 
allow to treat the different gauge fields as independent. 



dtd^r fi^iWI^ + —B") + A°(I^° 
mc 



2mc^ 



where S'[A^,W;,^/'] = J dt d^r C{Af„W;^,i:). For the 
usual electric and magnetic fields, we obtain the follow- 
ing equations of motion (Maxwell equations) via variation 
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w.r.t the Abelian gauge field components: 



Pmat. 



2mc- 



■dkeaMK = 0, (11) 



mc 



-^aikdkfJ'"' 



+ W^*^-'^ = 0' (12) 



with p^at. = e,p^ and J, = =^{^\d,^) - {d,ib^)^) - 
^Aitjj'^tjj. In equation (|lip. we can define a polarisation 
in terms of the Lagrange multipliers, = -f;^£akiK- 
This dielectric polarization is associated to bound charges 
and there appears then a contribution —dkPk to the to- 
tal charge density. This polarization also contributes the 
second equation of motion (|12|) through a term added 
to the ordinary current density, dtPi, as well as an ad- 
ditional term Cikadk'm"', with m° = describ- 
ing "Amperian" currents. The electric charge and cur- 
rent density are thus modified by the SO-|-Zeeman terms 
and the usual Maxwell equations become |28]: V • E = 
^^-V-P, and VxB-£oAio9tE = /^o[J + i9tP + V x m]. 

The variation with respect to the non-Abelian gauge 
fields leads to 

-Pmat. + ^^abcWlEl - -Ok^l + ^i^ = 0, (13) 

2c c 



hce,,kd,Bl + -dt^~X1^Q, (14) 



where instead of the field strength tensor GJ^^, we have 
defined the non-Abelian "electric" and "magnetic" fields 
(written in sanserif) 



B*^ = V X + ieabcW''"'x 



(15) 



Note that we have defined the matter spin density and 
superfluid spin current as 



Pmat. 



2m 



(V'V''(a,v) - (a.^^)T», 



(16) 
(17) 



where, as already discussed in ref. [5] , the presence of the 
gauge symmetry breaking term in the Lagrangian den- 
sity exactly compensates the "diacolor" contribution to 
the matter current and renders zero spin conductivity 
[24H26] , as expected from one electron Hamiltonians with 
linear SOI. We also recover the appearence of the "radia- 
tion" spin density and the "radiation" spin current density 
which would follow from the Noether theorem. 



Jt = 



ce,,kW^B%), 



(18) 
(19) 



in terms of which the corresponding Yang-Mills-Maxwell 
equations become 



V 



Jl VPmat. 



+ P?ad.-AO, 
1 



(20) 



V X B"" ~ ^dtE" = —(r + J" - X"), (21) 
he 

where the calligraphic characters denote the radiation con- 
tributions while roman fonts are used for the matter quan- 
tities. The occurrence of the radiation contribution to the 
conserved Noether current are directly linked to the exis- 
tence of the SO and Zeeman interactions and these terms 
are mandatory in order to satisfy the conservation of the 
total angular momentum. The matter contributions to 
this conserved current is the usual spin current associated 
to the free electrons, while the radiation part corresponds 
to the angular momentum carried by the lattice. 

In order to illustrate how the above contributions are 
at work in a simple case, let us again contemplate the 
case of a 2D electron gas with only uniform Rashba SO 
interaction, immersed in a perpendicular magnetic field. 
Using the notations k — 2ma/h and = eB/m, we 
write the components of the non-Abelian gauge field Q 
as = — Ky, = Kx and Wq — fl/c. The com- 
ponents of the non-Abelian electric and magnetic fields 
follow, and allow to compute the non vanishing compo- 
nents of the radiation spin density p] 



,3 

rad. 



-Tp-f^K and 

2c 



2c 



ad. 



spin current jT^ad 

— K^c^)x. The radiation spin current induces an 
electric polarization given by 



2mc3 ^2 



2 2\ 
K C )Z 



(22) 



Note that an expression similar to equation (j22p was ob- 
tained in Ref. j28| . The spin polarization generates a mag- 
netic moment whose radiation contribution is given by 



M 



eh 
2mc 



(23) 



The interesting result here is that a pure Zeeman term 
would only produce a spin precession, hence would not po- 
larize the electrons, but in association with SO interaction 
which acts with a perpendicular torque (relaxation of the 
spin polarization along the precession axis was first writ- 
ten phenomenologically by Landau and Lifshitz [29ll30] ). 
the spin orientation relaxes and a net polarization occurs. 

The electric polarization and the magnetic moment are 
physical observable quantities, and a gauge transforma- 
tion, if contemplated, would change the orientation of 
these fields, and hence the physical situation. 

From more general arguments, the necessity of a fixed 
gauge is clear when one notices that neither the matter 
(|16|17P nor radiation ()18|19|) spin polarizations or cur- 
rents are gauge invariant. In particular we have covari- 
ance of the matter current Jf — > Jf — eabcOt^Jf and an 
even more complicated transformation law for the radia- 
tion current [14] . Concerning the radiation contribution. 
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this is a general feature of non- Abelian gauge fields P"51IT^ . 
The case of the matter spin current is more perplexing 
and the present situation is due to the cancellation of the 
diacolor contribution. In the U{1) case, both the param- 
agnetic and the diamagnetic current densities are gauge 
dependent, but their sum is not. In the case of supercon- 
ductivity, where only one of the contributions survives (the 
diamagnetic one) , the current is made physical by the fact 
that the gauge is fixed [31]. We have in a sense a parallel 
situation here, where again, gauge fixing is required. 

Now the role of the Lagrange multipliers needs to be 
analyzed more completely. From the Yang-Mills-Maxwell 
equations, one can form a continuity equation describing 
the conservation of the total angular momentum density. 
Taking the divergence of Eq. (HU, and using Eq. ^IU\\ . one 
gets 

5t(pr„at. + P?ad.) + V ■ (J" + - at/i"^ + V • A^ (24) 

The Lagrange multipliers formally appear in the continu- 
ity equation, but physically, the l.h.s. already contains 
all possible sources of angular momentum in the problem, 
i.e. i) matter spin density (or all free electrons contri- 
butions from the conduction band encoded in the spinor 
ip)^ and ii) radiation contributions (or here angular mo- 
mentum transfercd to the lattice via SO interaction). The 
firts contribution accounts for conducting electrons contri- 
bution, since these latter electrons (e.g. from an s-band) 
do not contribute through an orbital angular momentum 
while the second contribution includes all external fields 
terms from the Lagrangian density ([3]). In the general 
case, there could also exist an angular momentum density 
associated to the ordinary U{1) fields [32], but since we 
did not incorporate any couplings of these fields to matter 
to allow for transfer of angular momentum, they should 
not enter the continuity equation. 

There could also exist angular momentum density as- 
sociated to the ordinary U{1) fields [32], but since we did 
not incorporate any couplings of these to matter to allow 
for transfer of angular momentum, they should not enter 
the continuity equation. So eventually the conservation of 
angular momentum in the system reads as 

5*(Pmat. + PU) + V • (J" + J^) = 0. (25) 

This equation is often discussed in the literature in the 
form 9tPmat + V • J° — torque density^ considering that 
we can arbitrarily modify the content of the matter spin 
current and that of the torque density in such a way that 
the equation is always satisfied. Our claim here is that 
such a freedom does not exist, and we can indeed inter- 
pret the radiation contributions, when written at the r.h.s 
of equation p5|. as a torque density, but the respective 
definitions of matter and radiation contributions are fixed. 

Physical discussion and measurements. A key 

point of our argument in this letter is the fact that due to 
the SU{2) symmetry of the gauge group for the SO inter- 
action, gauge transformations if performed, would change 



measurable physical quantities (as a result, these transfor- 
mations are of course not allowed). It is thus interesting 
to discuss in more detail possible measurements of such 
quantities. 

One of the primary quantities touched by such gauge 
transformations would be the spin polarization and spin 
current density. The electric polarization induced by the 
presence of spin current is an interesting feature for the de- 
tection of spin currents. In the same vein, the detection of 
charge imbalance was proposed a few years ago by Valen- 
zuela and Tinkham j33j . when they observed that due to 
the Onsager reciprocal relations between spin and charge 
currents, a measurable charge asymmetry is expected in 
the presence of a spin current. Let us also mention that 
the cubic dependence on the Rashba parameter of the di- 
electric polarization ([^^ induced by the radiation current 
should help in discriminating between matter and radia- 
tion contributions to the spin current density. 

An even more transparent example is given by the fa- 
mous experiments of the beginning of the XXth century 
based on gyromagnetic phenomena |341l35j . and among 
them the rotation of the sample in the Einstein-de Haas 
experiment. As we have argued, the continuity equation 
can be re-interpreted in terms of a torque density when all 
radiation contributions are shifted on the r.h.s. of the con- 
servation equation. This torque density resulting from a 
re- writing of the radiation polarization and current (in the 
case of the Rashba material treated as an example in this 
letter it would just be —dtp'^^^^ ~ j^Kil), it would thus be 
affected by a tentative gauge transformation. The "real- 
ity" of this torque density is demonstrated in the Einstein- 
de Haas experiment, where a magnetic field applied to an 
initially non magnetized ferromagnct induces a spin polar- 
ization, hence a spin angular momentum then transferred 
to the lattice to ensure conservation of the total angular 
momentum. As a result, one obtains an orbital motion of 
the lattice. 

A very illuminating mechanical analog of angular mo- 
mentum currents, relevant to the previous argument, is 
discussed in reference [T] : If one considers a solid cylinder 
made of any particular material and twists it by a net angle 
8 at one end, we generate a torsion. If one then attaches 
the two ends of the cylinder, thereby making a doughnut, 
the resulting torsion will be built into the material and 
cannot be relaxed except by plastic deformations. This 
is an example of a persistent angular momentum current 
and a source torque proportional to the gradient of the 
twisting angle. In the sense of the discussion before, once 
the material has been chosen and the end to end twisting 
angle is fixed then one has no freedom in changing the 
torque or the matter angular momentum current. This a 
clear example of a mechanical fixed gauge situation. 

Conclusions. — As we see, the variational definition 
of the currents in terms of derivatives of the Lagrangian is 
an unambiguous way to define a current, since its canoni- 
cal form does not depend on an ad hoc proposition. Nev- 
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ertheless, the variational principle by itself is not enough 
to solve the ambiguity of the spin current definition. To be 
consistent, we need also a symmetry broken formulation 
(or fixed gauge), otherwise, gauge transformations would 
change the matter (spin) and radiation (torque) content 
i.e. the matter currents are gauge covariant. A similar 
gauge symmetric scenario to the one presented here is 
materialized in superfluid condensates of neutral Bosons 
in the Helium 3 B phase, where one encounters a SU{2) 
Anderson- Higgs mechanism [5]. In reference [11], Tokatly 
argues on the fact that the variational definition of the 
current density leaves no room for any ambiguity on the 
definition of the spin current density. We would like to 
stress that this statement is only valid in the formulation 
discussed here, in a gauge symmetry broken formulation. 
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